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Mixing at Belle 

L. M. Zhang (on behalf of the Belle Collaboration) 
Syracuse University, Syracuse, New Yorl< 13244, USA 



We report the recent two results of D'^-D'^ mixing studies at Belle in — > K~^K /n~^n and — » Kgn~^n 
decays. The former measures the relative difference of the lifetimes ycp, giving the evidence of D^-D^ mixing; 
the latter measures the mixing parameters x and y. 
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1. Introduction 

Mixing phenomenon, i.e. the osciUation of a neutral 
meson into its corresponding anti-meson as a function 
of time, has been observed in the , , and most 
recently systems. This process is also possible in 
the D-meson system, but has not previously been ob- 
served. 

Mixing in heavy flavor systems such as that of B^ 
and B^ is governed by the short-distance box dia- 
gram. However, in the system this diagram is both 
GIM-suppressed and doubly-Cabibbo-suppressed rel- 
ative to the amplitude dominating the decay width, 
and thus the short-distance rate is very small. Conse- 
quently, mixing is expected to be dominated 
by long-distance processes that are difhcult to calcu- 
late; theoretical estimates for the mixing parameters 
X— (toi— m2)/r and y = (ri—r2)/2r range over two- 
three orders of magnitude [1]. Here, mi,m2 (ri,r2) 
are the masses (decay widths) of the mass eigenstates 
1^1,2) = p\D°) ± q\W), and r = (Fi + T2)/2. The 
parameters p and q are complex coefficients satisfying 
\p? + \q? = l. 

The general experimental method identifies the fla- 
vor of the neutral D meson when produced by recon- 
structing the decay D*+ -> D°7r+ or D*- D^n- 
[2] ; the charge of the accompanying pion identifies the 
D flavor. Because the energy release in D* decays is 
only ~ 6 MeV, the background is largely suppressed. 
The I?" decay time (t) is calculated via {l/p) x mj^o, 
where I is the distance between the D* and 1?° de- 
cay vertices and p is the D° momentum. The D* 
vertex position is taken to be the intersection of the 
D'^ momentum with the beamspot profile. To reject 
Z?^*' decays originating from B decays, one requires 
Pd* > 2.5 GeV, which is the kinematic endpoint. 



2. CP-eigenstates K+K~ and tt+tt" 



We have studied the decays to CP eigenstates 13° 
K'^K~ and Z3°—> tt+tt"; treating the decay-time dis- 
tributions as exponential, we measured the quantity 



where tj^-^+ and tj^+k- ^-re the lifetimes of D*^ — > 
7^-TT+ and D°^K+K-{or D^^tt+t:-) decays. It 
can be shown that ycp = ycos0 — ^AmxsiikP [3], 
where Am parameterizes CPV in mixing and </> is a 
weak phase. If CP is conserved. Am = </> = and 
ycp = y- This method has been used by numerous 
experiments to constrain y^p [4]. Our measurement, 
based on 540 fb~^ data, yields a nonzero value of ycp 
with > 3(7 significance [5]. We also searched for CPV 
by measuring the quantity 



At 
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this observable equals Ar = ■^AmV cos (j) — a; sine/) [3]. 

We reconstruct D*^ D^^^t decays and — > 
K^K~ , K^TT^, and tt+tt^. Candidate D° mesons are 
selected using two kinematic observables: the invari- 
ant mass of the decay products, M, and the energy 
release in the D*+ decay, Q = {Md* - M — m^+)c^ . 
According to Monte Carlo (MC) simulated distribu- 
tions of t, M and Q, background events fall into four 
categories: (1) combinatorial, with zero apparent life- 
time; (2) true mesons combined with random slow 
pions (this has the same apparent lifetime as the sig- 
nal) (3) decays to three or more particles, and (4) 
other charm hadron decays. The apparent lifetime of 
the latter two categories is 10-30% larger than tjjo. 

For the lifetime measurements, we select the events 
satisfying |AAf I/ctm < 2.3, \Q - 5.9 MeV| < 0.80 
MeV and at < 370 fs, where AM = M — rrijjo, and 
at is the decay time uncertainties calculated event- 
by-event. The invariant mass resolution ctm varies 
from 5.5-6.8 MeV/c^, depending on the decay chan- 
nel. The selection criteria are chosen to minimize the 
expected statistical error on ycp using the MC. We 
find 111 X 10^ K+K-, 1.22 x 10^ K~tt+ and 49 x 10^ 
TT+TT~ signal events, with purities of 98%, 99% and 
92% respectively. 

The relative lifetime difference ycp is determined 
by performing a simultaneous binned maximum like- 
lihood fit to the L»" ^ K+K-, D° K+tj- , D° 
TT+TT~decay time distributions. Each distribution is 
assumed to be a sum of signal and background con- 
tributions, with the signal contribution being a con- 
volution of an exponential and a detector resolution 
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function, 

dN/dt ■ 



'r ■ R{t~t')dt' + B{t). (3) 



The resolution function R(t — t') is eonstrueted from 
the normaUzed distribution of the decay time uncer- 
tainties at- The of a reconstructed event ideally 
represents an uncertainty with a Gaussian probabil- 
ity density: in this case, bin i in the at distribution 
is taken to correspond to a Gaussian resolution term 
of width (Ti, with a weight given by the fraction fi 
of events in that bin. However, the distribution of 
"pulls", i.e. the normalized residuals (tree — igen)/^* 
(where tree and tgen are reconstructed and generated 
decay times), is not well-described by a Gaussian. We 
found that this distribution can be fitted with a sum 
of three Gaussians of different widths cr^"'' and frac- 
tions Wk, constrained to the same mean. Therefore, 
we choose the parameterization 



= l^/iI^«'fcG(t-t';aifc,io), (4) 



i=l 



i=l 



with aik — Sfecr^^'^CTi, where the Sk are three scale fac- 
tors introduced to account for differences between the 
simulated and real crj^"^\ and to allows for a (common) 
offset of the Gaussian terms from zero. 

The backgroimd B{t) is parameterized assuming 
two lifetime components: an exponential and a 6 func- 
tion, each convolved with corresponding resolution 
fimctions as parameterized by Eq. (4). Separate B{t) 
parameters for each final state are determined by fits 
to the t distributions of events in M sidebands. The 
MC is used to select the sideband region that best re- 
produces the timing distribution of background events 
in the signal region. 

Fitting the K~tt~^, K~^K~, and tt+tt" decay time 
distributions (Figs. l(a)-(c)) shows a statistically sig- 
nificant difference between the K~'it~^ and h^h~ life- 
times. The effect is visible in Fig. Id, which plots 
the ratio of event yields A^^+^./AT^^ as a function of 
decay time. The fitted lifetime of D'^ meson in the 
K^n^ final states is 408.7 ± 0.6 fs, which is consis- 
tent with the PDG value [6] (and actually has greater 
statistical precision). We measure 



Vcp = (1.31 ± 0.32 ± 0.25)% , 



which deviates from zero by 3.2(T. The systematic 
error is dominated by uncertainty in the background 
decay time distribution, variation of selection criteria, 
and the assumption that is equal for all three final 
states. The analysis also measures 



(0.01 ±0.30 ±0.15)%, 



(6) 



which is consistent with zero (no CPV). The sources 
of systematic error for Aj. are similar to those for y^p . 




Figure 1: Projections of the decay-time fit superimposed 
on the data for D° -^K+K', D° ^ K--K+, and D° 
'K^'K~ decays. The hatched area represents the back- 
ground contribution. Plot (d) shows the ratio of decay- 
time distributions for {K+K~ + tt+tt") and D° 
K~n^; the solid line is a fit to the points. 



3. Dalitz Plot Analysis of 



The time dependence of the Dalitz plot for D° 
Kg TT+TT^ decays is sensitive to mixing parameters x 
and y without ambiguity due to strong phases. For 
a particular point in the Dalitz plot (to+, rri^), where 
m_,_ = m{KgTT~^) and m_ = m{Kg Tr~), the overall 
decay amplitude is 



where e^^2)it) = e"(""i.2+^i'2/2)t. The first term 
represents the (time-dependent) amplitude for 
KgiT~^Tr^, and the second term represents the ampli- 
tude for D^^D^^K'^TT+n-. Taking the modulus 
squared of Eq. (7) gives the decay rate or, cquiva- 
(5) lently, the density of points p(m'^,m^;t). The result 



contains terms proportional to cosh(yrt), cos(xri), 
and sin(a;rt), and thus fitting the time-dependence of 



/,t) determines x and y. This method was 
developed by CLEO [7]. 

To use Eq. (7) requires choosing a model for the 

decay amplitudes A^^^ —g{m'^,m'^). This is usually 

taken to be the "isobar model" [8], and thus, in ad- 
dition to X and y, one also fits for the magnitudes 

and phases of various intermediate states. Specifi- 
cally, Ajjo{ml,ml) = T.jOije'^^A 



where 5^ is a 
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strong phase, Aj is the product of a relativistic Breit- 
Wigner function and Blatt-Weiskopf form factors, and 
the parameter j runs over all intermediate states. This 
sum includes possible scalar resonances and, typically, 
a constant non-resonant term. For no direct CPV, 
A—g{m'^,mt) — .4^0 (jTi?., m^); otherwise, one must 
consider separate decay parameters (a^ , S^) for D*^ de- 
cays and {aj,Sj) for D° decays. 

We have fit a large D'^ tt+tt" sample selected 

from 540 fb^^ of data [9]. The analysis proceeds in 
two steps. First, signal and background yields are 
determined from a two-dimensional fit to variables 
M{Kmr) and Q = MiTT^Kmr) - M{K7nr) - m^+. 
Within a signal region \M {Knn) — m^o | < 15 MeV/c^ 
and \Q — 5.9 MeV| < 1.0 MeV (corresponding to in 
resolution), there are 534000 signal candidates with 
95% purity. These events are fit for x and y; the 
(unbinned ML) fit variables are m^, m^, and the 
decay time t. Most of the background is combina- 
toric, i.e., the Z?" candidate results from a random 
combination of tracks. The decay-time distribution 
of this background is modeled as the sum of a delta 
function and an exponential function convolved with 
a Gaussian resolution function, and all parameters 
are determined from fitting events in the sideband 
30 MeV/c2 < iMiKiTTT) - m^o| < 55 MeV/c^. 

The results from two separate fits are listed in Ta- 
ble I. In the first fit CP conservation is assumed, i.e., 
q/p—1 and Aj^g{m'^,m^) = ^^o('Ti?_, m^). The free 
parameters are x,y, r^o , some timing resolution func- 
tion parameters, and decay model parameters (a^, S^.). 
The results for the latter are listed in Table II. The 
results for x and y indicate that x is positive, about 2a 
from zero. Projections of the fit are shown in Fig. 2. 
The fit also yields r^, = (409.9 ± 1.0) fs, which is 
consistent with the PDG value [6] (and actually has 
greater statistical precision). 

Table I Fit results and 95% C.L. intervals for x and y, 
from analysis of — > Kg tt+tt" decays. The errors are 
statistical, experimental systematic, and decay-model sys- 
tematic, respectively. 



Fit Par am. 


Result 


95% C.L. inter. 


No X (%) 
CPV y{%) 


0.80 ± 0.29 in? iow 
0.33 ± 0.24 inilSfs 


(0.0,1.6) 
(-0.34,0.96) 


CPV X (%) 

y(%) 
\<i/p\ 


0.81 ±0.30 lg:J?lgfe 

0.37 ±0.25 tViltVol 
0.86 «;^°iroU 0.08 

1/1 +16+5+2 
-18 -3 -4 


|a;| < 1.6 

\y\ < 1-04 



For the second fit, CPV is allowed and the 13° 
and D samples are considered separately. This in- 
troduces additional parameters \q/p\, ATg{q/p) = cf), 
and {aj,Sj). The fit gives two equivalent solutions, 
{x,y,(p) and (— x, — y, (/> + tt). Aside from this pos- 



Table II Fit results for iCg tt+tt" Dalitz 

plot parameters. The errors are statistical 

only. The fit fraction is defined as the ratio 
of the integral J la^ Ar{m^ ,m^)\^ drrt^ dni\. to 
/ I Er=i "-r e'^^^Arim^. , m\)f dml dm\. 



Resonance 


Amplitude 


Phase (deg) 


Fit fraction 


if* (892) ~ 


1.629 ± 0.006 


134.3 ± 0.3 


0.6227 


Kq (1430) 


1 _i_ n AO 


A n _i_ A 

— u.y ± u.o 


A A'TO /I 

0.0/24 


K2 (1430) 


0.87 ± 0.02 


—47.3 ± 1.2 


0.0133 


K (1410) 


0.65 ± 0.03 


111 ± 4 


0.0048 


K (1d80) 


O.dO ± 0.25 


147 ± 29 


A AAAO 

0.0002 


A (892)^ 


0.152 ± 0.003 


—37.5 ± 1.3 


0.0054 


(1430)^ 


0.541 ± 0.019 


91.8 ± 2.1 


0.0047 


A:|(1430) + 


0.276 ± 0.013 


-106 ± 3 


0.0013 


K*{IA:10)^ 


0.33 ± 0.02 


— 102 ± 4 


0.0013 


it"* (1680)+ 


0.73 ±0.16 


103 ± 11 


0.0004 


p(770) 


1 (fixed) 


(fixed) 


0.2111 


a;(782) 


0.0380 ± 0.0007 


115.1 ± 1.1 


0.0063 


/o(980) 


0.380 ± 0.004 


-147.1 ± 1.1 


0.0452 


/o(1370) 


1.46 ±0.05 


98.6 ±1.8 


0.0162 


/2(1270) 


1.43 ±0.02 


-13.6 ± 1.2 


0.0180 


/9(1450) 


0.72 ±0.04 


41 ±7 


0.0024 


CTl 


1.39 ±0.02 


-146.6 ±0.9 


0.0914 


0-2 


0.267 ±0.013 


-157 ±3 


0.0088 


NR 


2.36 ±0.07 


155 ± 2 


0.0615 




Figure 2: Projection of the unbinned ML fit superimposed 
on the data for D^^iCS TT+TT decays. In (d), the hatched 
area represents the combinatorial background contribu- 
tion, and the lower plot shows the ratio of decay-time 
distributions for events in the it* (892)+ and A" (892)" 
regions, where sensitivity to {x, y) is highest. 
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siblc sign change, the effect upon x and y is small, 
and the results for \q/p\ and (p are consistent with 
no CPV. The sets of Dalitz parameters {a^,S^) and 
(a,., (5,.) are consistent with each other, indicating no 



direct CPV. Taking a, 



: Qj and Sj - 



(i.e., no direct 



0.95 



+0.22 
-0.20 



CPV) and repeating the fit gives \q/p\ 
and^= (-2+}°)°. 

The dominant systematic errors are from the time 
dependence of the Dalitz plot background, and the 
effect of the pjj, momentum cut used to reject D*'s 
originating from B decays. The default fit includes tttt 
scalar resonances ai and a2 ; when evaluating system- 
atic errors, the fit is repeated without any tttt scalar 
resonances using if- matrix formalism [10]. The influ- 
ence upon X and y is small and included as a system- 
atic error. 

The 95% C.L. contour for (x, y) is plotted in Fig. 3. 
The contour is obtained from the locus of points where 
— 21n£ rises by 5.99 units from the minimum value; 
the distance of the points from the origin is subse- 
quently rescaled to include systematic uncertainty. 
We note that for the CPT^-allowed case, the reflec- 
tions of the contours through the origin are also al- 
lowed regions. 
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Figure 3: 95% C.L. contours for (x,y): dotted (solid) 
is statistical (statistical plus systematic) contour for no 
CPV\ dashed-dotted (dashed) is statistical (statistical 
plus systematic) contour allowing for CPV . The point 
is the best-fit value for no CPV . 
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